Theoretical study of the accuracy limits for the optical resonance frequency measurements. 
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The principal limits for the accuracy of the resonance frequency measurements set by the asymmetry of 
the natural resonance line shape are studied and applied to the recent accurate frequency measurements in the 
two-photon ls-2s resonance and in the one-photon ls-2p resonance in hydrogen atom. This limit for ls-2s 
resonance is found to be ~ 10 Hz compared to the accuracy achieved in experiment ±46 Hz. In case of 
deuterium atom the limit is essentially larger: 10 -2 Hz. For ls-2p resonance the accuracy limit is 0.17 MHz 
while the uncertainty of the recent frequency measurement is about 6 MHz. 



The resonance approximation which reduces the descrip- 
tion of the spectral line profile to the Lorentz contour is the 
cornerstone of all the modern transition frequency measure- 
ments in the resonance experiments. This approximation is 
employed for the extraction of the resonance frequency values 
from the experimental data. Here we address only the natural 
line profile, disregarding Doppler and collisional broadening. 
We assume that the latter effects can be minimized below the 
level of the natural line width. In the resonance approxima- 
tion the line profile is described by two parameters, resonance 
frequency uf es and the width V and is symmetric with respect 
to ui res . In the early paper by F. Low 0JJ] it was pointed out 
that the resonance approximation is valid only up to a certain 
limit of accuracy which is defined by nonresonant (NR) cor- 
rections. Beyond this limit the line profile acquires asymmet- 
ric form and the definition of u> res becomes nonunique. More- 
over, the NR corrections are, in principle, process-dependent, 
so the frequency value also will differ for different excitation 
processes. This means that the measurement of frequency is a 
well defined procedure only up to some limit, which we define 
in our paper. Going beyond this limit, it should be replaced by 
the direct measurement of the spectral line profile, which then 
becomes the only quantum-mechanical observable. 

During the last few years the NR corrections were recon- 
sidered in a series of theoretical works |2J]-[l5[]. This reconsid- 
eration was triggered by the new extra accurate optical reso- 
nance experiments In the experiment [6] the natural 
line profile for the Lyman-a ls-2p one-photon resonance in 
hydrogen was observed for the first time. In Q7TJ , the most 
accurate measurements in the optical region were performed 
for the two-photon ls-2s resonance in hydrogen. In the lat- 
ter experiments the absolute accuracy of the frequency mea- 
surements was as high as ±46 Hz, or 10~ 14 in relative units. 
However, the scheme of these experiments causes serious dif- 
ficulties for the theoretical description of NR corrections and 
hence for the definition of the accuracy limits in this case. In 
[0, A, the H atoms were excited from the ground Is state to 
the 2s state via two-photon absorption in the space region free 
of external fields. After this the excited atoms moved during 
approximately 10~ 3 s to another space region where they meet 
a weak electric field. In this field 2s and 2p states are mixed 



and the atoms decay via the ordinary 2p-ls transition. This ra- 
diation is detected and provides the necessary information for 
the extraction of ls-2s frequency value from the experimental 
data. In the rest frame of an atom it looks like the excitation 
occurs in the absence of the electric field and then the electric 
field is turned on (delayed decay). In terms of QED the ini- 
tial and final states of an atom are described then by different 
Hamiltonians Hi n and H out . 

The QED theory with different in- and out-Hamiltonians 
was developed by Fradkin, Gitman and Shvartsman (see book 
|0]). This theory we will use below for the description of NR 
corrections in the two-photon ls-2s resonance experiment. 
The earlier attempts 13D-I1J] to describe the NR corrections 
in these experiments with standard QED methods cannot be 
considered as reliable. 

The major contribution to NR corrections arises from the 
interference between the resonant and nonresonant terms in 
the expression for the scattering amplitude 10. In 10 only the 
nonresonant terms with the same symmetry as the resonant 
one were included. As it was pointed out in [0] the most im- 
portant contribution arises from the interference between the 
resonant term and nonresonant terms with another symmetry: 
the NR contribution of 2p 3 / 2 state to Is — 2pin resonance 
was considered as an example. It was argued that this contri- 
bution can survive in differential (with respect to the angles) 
cross-section for the resonance photon scattering on an atom. 
In the present work we will show that the NR contributions 
of this type survive both in differential and total cross-section 
provided that the resonance line shape is fully natural. More- 
over, in many cases the dominant NR contributions which de- 
fine the accuracy limits for the frequency measurements arise 
from the neighbor hyperfine (HF) components. These com- 
ponents have different symmetry (values of the total angular 
momentum F of an atom) compared to the HF component 
which is used in the basic transition (see Fig. 1). The contri- 
bution of HF neighbor components was missing in the earlier 
investigations 10-10]. 

We start with the simplest case of the elastic photon scatter- 
ing on the hydrogen atom in the absence of the electric field. 
The amplitude for this process can be written in a standard 
formlHoll: 
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where the quantum numbers njlm represent the standard 
set of one-electron quantum numbers for the electron in the 
hydrogen atom. The photons are characterized by the total 
angular momentum J, its projection M and the orbital an- 
gular momentum L. The latter defines parity, or the type of 
the photon: electric or magnetic. The operator Ajlm corre- 
sponds to photon absorption, A* JLM denotes photon emission. 
In Eq. (Q} the additional term where the operators A* and 
A are interchanged and the sign of uj is reversed, is omitted. 
This term vanishes in the resonance approximation and does 
not contribute to the leading NR corrections. The resonance 
condition for the frequency is uj = E ni j 1 i 1 — E no j i and in 
the resonance approximation in Eq. (Q]i only the terms with 
njl = nxjili are retained. 

To avoid the singularity in Eq. (Q]l we insert in the de- 
nominator the width of the excited state T ni j 1 i 1 . The regu- 
lar method for this insertion within the framework of QED is 
described in H, d, d, 10. 



With this description of the resonant process the NR correc- 
tions arise as the other terms of the expansion in Eq. (UJ with 
njl 7^ nxjih- Performing the summations over all the angu- 
lar momentum projections, we obtain for the process proba- 
bility in the resonance approximation: 
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where W ab , W em are the standard absorption and emis- 
sion probabilities per time unit and S 1 (joJij2"0 = 
( jo ji J } 

i J J a >■ Taking into account the NR correction due 
a { 32 jo J J 

to the closest to n\j\l\ neighbor level n 2 j2h results in 
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Here A em , A ab are the reduced emission and absorption 
amplitudes (reduced matrix elements of the photon emission 
and absorption operators). These matrix elements do not 
depend on the angular momentum projections. The factor 
S ll {jojij2J) does not vanish for j% ^ j\, so that the dom- 
inant NR corrections due to the interference of the resonant 
and nonresonant terms with different symmetry 7^ jl) 
can exist. They exist both in differential (with respect to the 
angles) and total probability; Eq. (01 describes the total proba- 
bility. The differential one would depend on the angle (rfj, n}) 
where ni, rif are the velocity directions for the absorbed and 
emitted photons. The vector ni is fixed by the incident laser 
beam and rif by the position of the detector. 

Now we turn to the case of the ls-2s two-photon absorption 
and delayed decay in electric field. The formalism of QED 
with different in- and out-Hamiltonians developed in |9|] pur- 
sues the more complicated task: the creation of the electron- 
positron pairs in the strong electric field. We will apply this 
formalism to the case of a weak electric field. However, the 
situation in our case is also non-perturbative: due to the ad- 
mixture of the 2p state into the 2s state the emission probabil- 
ity changes by 8 orders of magnitude. The change of the other 
atomic characteristics (Stark shift of the energy levels, Stark 
splitting) we will consider as negligible. The criterion of the 
weak field will be e < e c = 475 V/cm, where e is the strength 



of the electric field. In the field e = e c the 2s and 2p levels are 
100% mixed Q. 

The Fradkin-Gitman-Shvartsman (FGS) theory y9[] operates 
with two complete sets of eigenfunctions belonging to in- 
and out-Hamiltonians. This theory follows the standard QED 
approach in generalized form: S-matrix, field operators in 
the Fock space, 4-dimensional perturbation expansion for the 
S-matrix elements, Wick theorem and Feynman graph tech- 
niques. Actually the unique new element that we will have to 
use is the generalized FGS electron propagator. This propa- 
gator connects two vertices, which are described by in- and 
out-Hamiltonians, respectively. 

In our case this means the absence or presence of the elec- 
tric field. The FGS propagator looks like 

s FGS {x lX2 ) = e(t! - 1 2 ) V^iK^nO^w 

fi,n(-E«;ra>0) 
n,n(B fi; „<0) 

where ipa(x) are the solutions of the Dirac equation for the 
electron in the field of the nucleus and the external electric 
field, i/j n (x) are the solutions with zero external field; En, 
E n are the corresponding eigenvalues. The matrix ui nn in the 
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weak-field limit reduces to an overlap integral 

uJnn = J ip£ (x)ip n (x)dx =< n\n > (5) 

In the nonrelativistic limit, evidently valid for the neutral 
hydrogen atom, we replace the Dirac wave functions by 
Schrodinger ones and omit the negative-energy contribution 
in Eq. ©. 

Within the FGS theory, the probability of absorption of two 
equivalent laser photons with frequency ui' by the electron in 
the hydrogen atom in its ground state with the subsequent de- 



layed decay in external electric field in the resonance approx- 
imation (n = a, n = a') looks like 
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where W^ ah21 ^ is the two-photon absorption probability, S a 2> 
is the angular factor similar to 5 l7 in Eq. (0, and the reso- 
nance condition is E a > — E a — 2uj' . 

Taking into account the NR correction due to the closest to 
a' neighbor level a" results in an expression similar to Eq. ([3]) 
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where A^ b a 2j \ A^ 2 ^ ^ tne reduced two-photon ab- 
sorption amplitudes. 

For moving further we have to choose the procedure for de- 
termination of the resonance photon frequency. In fl-lHtl the 
evaluation of the maximum value of the frequency distribu- 
tion was used for this purpose. As it was shown in yfl any 
other procedure (e.g. finding a "center of gravity" for the line 
profile) would give the result quite close to the choice formu- 
lated above. In case of the Lorentz profile all the methods of 
defining uf es give the same result Lo r a es = E a i — E a . With our 
choice the NR correction will look like 



Su> NR = u max 
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where the frequency value, corresponding to the max- 

imum of the frequency distribution. 

The first example, that we will consider, is the two-photon 
(F — 1) + 27 — > 2s 1 / 2 {F — 1) transition in hydro- 
. The value K = 1 for the total angular mo- 
mentum K of the two-photon system with equal photon fre- 
quencies is strictly forbidden II 1 Oil . This concerns exactly the 
case of the two-photon absorption of laser photons in experi- 
ment 10], |@] (see Fig. 1). According to this picture, the main 
NR contribution arises from the transition lsi/ 2 (F = 1) — » 
2pi/ 2 (F — !)■ F° r deriving the NR correction we have to use 
an expression (0 where we have to set a = \s\/ 2 (F = 1), 
a' = 2s 1 / 2 (F — 1) and a" = 2pi/ 2 (F — 1). In a weak elec- 
tric field Ipa ^ Ipa, ^a< — Vv + r]4>a" , Ipa" - Ipa" ~ T)1p a ' 

and the overlap integrals are < a\a >~< a"\a" >~ 1. 
Here 77 = \AEs\/AEl is the Stark shift to the Lamb shift 
ratio. This ratio is 77 = 1 for a field e = e c . For deriv- 
ing the NR correction we set in Eq. (0 E a > — E a = ^5 es , 
E a „ -E a = ujl es - AE L , where AE L = E a „ - E a , « 10 3 
MHz is the Lamb shift. The width T a > in case of the exper- 
iment [0], 10] is determined by the experimental set up (the 
time delay before the excited atoms enter the electric field re- 



wave functions ipa, rp~, and %/}-„ in the emission amplitudes 
yields: At™, « ?7(r 2p ) 1/2 , « (r 2p ) 1/2 , where r 2p is 

the width of 2p state. Then, evaluating the maximum value 
of Eq. (0 with respect to oj' and defining the NR correction 
according to Eq. (0 we find 
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where SI, (F, F') is the angular factor defined for the transi- 
tion between the two hyperfine sublevels. In the present case 
F = F' = F" = 1, and the corresponding angular fac- 
tors are equal: S*2' (F, F') = Sll„(F,F") = 11/18. The 
probability for the two-photon transition W^J^ 1 * s verv 
well known. An accurate nonrelativistic value for this transi- 



tion was obtained in (14(1: W 1 
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gion) and is equal to T a i = T e 



1 kHz. Insertion of the 



a.u. Here a is the fine structure constant, Z is the nuclear 
charge and it is assumed that aZ <C 1. The probability 
W\] 2pin was evaluated recently numerically for all Z val- 
ues (1 < Z < 100) lO and analytically for aZ < 1 
M- The result is W^ 2pi/2 = W%%^ + = 
2.907 • 10~ 5 (aZ) 8 + 3.69 • 10~ 6 (aZ) 8 (a.u.). 

The dependence r/^ 1 in Eq. (0 cannot be continued to zero 
field; the meaningful limit is set by the field where the de- 
cay rate of 2s level due to the admixture of 2p state becomes 
equal to the natural decay width of 2s state. This limiting field 
strength would be so small, that it cannot be used in a real 
experiment. Inserting all the numbers in Eq. (O and taking 
i] = 0.1 which corresponds to the weak electric field e — 47, 5 
V/cm, we obtain the final result \8lo nr \ « 10~ 5 Hz. This ac- 
curacy limit is still far from the recent inaccuracy estimate in 
the experiments H, 01 : ±46 Hz. However, it is important 
to notice that the situation in hydrogen seems to be rather for- 
tunate for the accurate resonance frequency measurement due 
to the absence of the transition to another (F — 0) hyperfine 
sublevel of 2s state (see Fig. 1). For comparison, in the deu- 
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terium where the total atom angular momentum values for the 
Is and 2s levels are F — 1/2, 3/2 respectively, this transition 
is allowed and NR correction by order of magnitude is 



F=1 



\Slo 



NR\ 



r 2 

exp 



l 

4'AEhfs 



(10) 



100 



where AE HFS = E 2si/2 ( F=3/2) ~ ^2s 1/2 (f=i/2) 
MHz is the hyperfine-structure interval for 2s level. Taking the 
same T exp value as in J3l, 0] we would have \Slu nr \ ~ 10~ 2 
Hz which is 3 orders of magnitude larger than for hydrogen. 
This is not so far from the accuracy limit of about 0. 1 Hz that 
was considered in [8] as achievable in the future with the use 
of colder hydrogen atoms. 

A quite different situation arises for the ls-2p resonant 
experiment in hydrogen Jg]. In (0] the resonance ls(F = 
1) — » 2p 3 / 2 {F — 2,1) was measured. The hyperfine 
structure for the 2p 3 / 2 level was not resolved since T 2p > 
AE HFS {2p 3 /2): T 2p « 100 MHz, AE HFS (2p 3/2 ) = 
E(2p 3/2 ,F = 2)~E{2p 3/2l F = 1) = 23.7 MHz. This 
is a typical case of overlapping resonances for two hyperfine 
sublevels. However, due to the presence of interference terms 
(see Eq. (0l) the line shape deviates from the overlap of two 
Lorentz profiles and can be presented by an expression: 
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where uj[ es = E(2p 3/2 ,F = 2) - E(ls 1/2 ,F = 1), u r 2 
E(2p 3 / 2 , F = 1) — E(lsi/ 2 , F= i), /, g are the angular factors, 
similar to S* 17 in Eq. ©. These angular factors play an impor- 
tant role in defining the line shape Eq. ( fTTb : /(1,2) : /(1,1) : 
g(l, 2, 1) = 181 : 1 : 0.307. Thus the line shape exhibits the 
one-peak structure (the larger peak fully screens the smaller 
one). The deviation from the overlap of two Lorentz profiles 
is of order 0.307/181 « 0.17%. Then the absolute accuracy 
limit for determination of ls(F = 1) — » 2p 3 / 2 transition fre- 
quency can be estimated as T 2p - 0.0017 « 0.17 MHz. The ex- 
istence of the interference terms that distort the closely lying 
resonances is well known (see for example, lfl7ll ). However, 
in (vi\\ this distortion was included in the error bars for the 
transition frequencies not distinguishing between other con- 
tributions of purely technical origin, such as the laser inten- 
sity distribution and the spread of atomic velocities. Nowa- 
days, when the two latter effects in practice can be diminished 
completely (e.g. by the use of colder atoms) it becomes more 
useful to introduce such quantity as an "absolute limit" for the 
frequency determination in the resonance experiments. Actu- 
ally any resonance experiment pretending on the utmost ac- 
curacy should contain fittings of the experimental data with a 
suitable theoretical line shape including the effects of interfer- 
ence, originating from NR corrections. In case of Is — 2p 3 / 2 
transitions Eq. (11) provides such a theoretical line shape, de- 
rived from first principles of QED. 
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Fig. 1. Scheme of the levels for the two-photon ls-2s transition. The vertical double 
lines denote the two-photon transitions. The K numbers denote the total angular 
momentum for a two-photon system, possible for different two-photon transitions. 
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